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EXTREMAL DICHOTOMY FOR UNIFORMLY HYPERBOLIC 

SYSTEMS 
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MARK HOLLAND, AND MATTHEW NICOL 


Abstract. We consider the extreme value theory of a hyperbolic toral automorphism 
T : T 2 —> T 2 showing that if a Holder observation <f> which is a function of a Euclidean- 
type distance to a non-periodic point C, is strictly maximized at £ then the corresponding 
time series {<j>o T 1 } exhibits extreme value statistics corresponding to an iid sequence of 
random variables with the same distribution function as <p and with extremal index one. 
If however <j> is strictly maximized at a periodic point q then the corresponding time-series 
exhibits extreme value statistics corresponding to an iid sequence of random variables with 
the same distribution function as <p but with extremal index not equal to one. We give a 
formula for the extremal index (which depends upon the metric used and the period of q). 
These results imply that return times are Poisson to small balls centered at non-periodic 
points and compound Poisson for small balls centered at periodic points. 


1. Introduction 

Suppose we have a time-series {X, } of real-valued random variables defined on a probability 
space (Q,P). Extreme value theory concerns the statistical behavior of the derived time 
series of maxima M n := max{Xi,..., X n }. There is a well-developed theory for iid random 
variables pn m nn and these references also consider the case of dependent random 
variables under certain assumptions. In the setting of ergodic maps T : X —> X of 
a probability space ( X , //) modeling deterministic physical phenomena we are concerned 
with time-series of the form {(po T 1 } where (p : X — > M is an observation of some regularity, 
usually Holder. In the literature, as in this paper, it is customary to assume that (p is 
a function of the distance d(x,p) to a distinguished point ( for some metric d so that 
(p(x) = f(d(x,Q) f° r x € A". For convenience the observation <p(x) = — logd(x,£) is often 
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used, but scaling relations translate extreme value results for one functional to another 
quite easily pm. Since the function £> o T l (x) detects the closer T l (x) is to £, there is a 
natural relation between extreme value statistics for the time series {c fioT *} and return time 
statistics to nested balls about £ (31 HO] - I n fact, if (X,d) is a Riemannian manifold and 
H is absolutely continuous with respect to the volume measure then return time statistics 
may be deduced from extreme value statistics and vice-versa, so in some sense they are 
equivalent though proved with different techniques. 

Due to the close relation with return time statistics one expects a difference in the extreme 
value theory of observations centered at periodic points to those centered at generic points. 
We refer to the early paper of Hirata [18] • What is perhaps surprising is that in many cases 
of uniformly expanding smooth one dimensional systems there is a strict dichotomy between 
the behavior of observations maximized at periodic points and other non-periodic points. 
Suppose (j) is maximized at a point £ e X and define a sequence of constants u n by the 
requirement that lim^^ n/i(0 > u n ) = r, where r G M. ft has been shown (7( [US 22j that 
for a variety of smooth one-dimensional chaotic systems this implies that lim r ,^ oc /i(M„ < 
«„) = e~ r}r where 0 < $ < 1 is called the extremal index and roughly measures the 
clustering of exceedances of the maxima (more precisely $ _1 is the expected number of 
observed exceedances if one exceedance is observed), ft turns out that d = 1 if £ is not 
periodic and •& < 1 otherwise (a precise expression for d may often be given as a function 
of the multiplier derivative of T along the periodic orbit). For higher dimensional smooth 
chaotic systems a similar dichotomy is expected (the presence of discontinuities complicates 
matters Cl) but so far has not been established. In this paper we establish such a dichotomy 
for hyperbolic toral automorphisms, give an expression for the extremal index and find that 
it also depends upon the metric used. This has immediate implications for the return time 
statistics. We have Poissonian return time statistics to non-periodic points and a compound 
Poissonian distribution to periodic points (the form of the compound distribution is related 
to the extremal index, as well). In fact, our results are given in terms of convergence of 
point processes of rare events to the standard Poisson process, in the non-periodic case, and 
to a compound Poisson process, in the periodic case, where the multiplicity distribution 
is also seen to depend on the metric used. We remark that the convergence of the point 
processes is actually stronger than establishing a Poissonian distributional limit for the 
number of visits to shrinking target sets. 

Earlier results for invertible dynamical systems include work by Denker, Gordin and 
Sharova, [5], who proved that, for toral automorphisms, the number of visits to shrinking 
neighbourhoods of non-periodic points converge (when properly normalised) to a Poisson 
distribution and work by Dolgopyat [6j, who also proved Poissonian return time statis¬ 
tics for non-periodic points in uniformly partially hyperbolic dynamical systems and noted 
possible generalizations to the case of periodic points. There have been a variety of re¬ 
cent results on extreme value and return time statistics for invertible dynamical systems, 
both uniformly and non-uniformly hyperbolic, but restricted to observations maximized at 
generic points 


EXTREMAL DICHOTOMY FOR UNIFORMLY HYPERBOLIC SYSTEMS 


3 


2. The setting 


Let L : M 2 —> M 2 be a 2 x 2 matrix with integer entries, with determinant equal to 1 and 
without eigenvalues of absolute value 1. The eigenvalues of L are irrational numbers A and 
1/A satisfying |A| > 1; the corresponding eigenspaces, E s and E u , are lines with irrational 
slope and the family of lines parallel to these subspaces is invariant under L. Moreover, 
L contracts vectors in E s and expands in E u . As L(Z 2 ) = Z 2 , this linear map induces an 
invertible smooth transformation T on the flat torus T 2 = M 2 /Z 2 through the canonical 
projection. The subspaces E s and E u project on dense curves which intersect densely and 
are T invariant; that is, the torus is the homoclinic class of the (unique) fixed point of T 
and a global attractor. 

After finding a norm in M 2 that makes L a uniform contraction in E s and a uniform 
expansion in E u , we project the Riemannian metric generated by this norm to T 2 and 
consider the invariant splitting at each point into one-dimensional affine spaces parallel to 
the subspaces E s and E u . This way, T becomes an Anosov diffeomorphism: the hyperbolic 
splitting at each point of T 2 is obtained by translating to that point the eigenspaces of 
the matrix. Through each point x G T 2 there is a local stable manifold W*(x) of length 
e and similarly a local unstable manifold l'T“(x). Besides, T is topologically mixing and 
its periodic points (all saddles) are dense, even though for each q G N the set of periodic 
points with period q is finite. 

As the determinant of L is equal to 1, the Riemannian structure induces a Lebesgue measure 
m on T 2 which is invariant by T. The probability measure m is Bernoulli and it is also 
the unique equilibrium state of the Holder map x G T 2 —» ip u (x) = — log | det DTru( x ) \, 
whose pressure is zero. This means that the metric entropy of m is given by h m (T) = 
— f <p u dm = log (A). 

Unless said otherwise, we assume that we are using the Euclidean metric in M 2 , which we 
project to T 2 and denote by d(-, •). More precisely, if (x,y) denote the local coordinates of 
the point z G T 2 then 

d{z, 0) = sjx 2 + y 2 . (2.1) 

Given a non-empty set A and a point x, we denote by A the closure of the set A and define 
the distance of re to A by dist(rc, A) := inf {d(x,y) : y G A}. 

Remark 2.1. Let A be a set homeomorphic to a ball whose boundary is a piecewise dif¬ 
ferentiable Jordan curve. For £ > 0, let D e := {x G T 2 : dist(rr,A) < e}. Then 
m{D e \ A) — e\dA\ + o(e), where | • | denotes the length of the curve dA. 


Recall that a continuous function ip : T 2 —y M is a-Holder if there are constants C > 0 
and a > 0 such that | (p(x) — <p(y)\ < Cd(x,y) a . For a G ]0,1], the space of a-Holder 
continuous functions with the norm 


IMI 


a 


IMIc° + sup 

x^y 


MX) - <p(y) I ) 


d{x,y) a 
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is a Banach space. And, given k G N and a G ]0,1[, any a-Holder continuous function 
p on the torus admits a Lipschitz approximation pi<. such that \\p — ^fc||oo — 0(k~ TZ ^). 
From the computations in (21 Section 1.26] one can deduce that the diffeomorphism T has 
exponential decay of correlations with respect to a-Holder functions, for any a G ]0,1]. 
More precisely, there exist constants C > 0 and 0 < 9 < 1 such that, for all ct-Holder maps 
(j) and ip, we have 


' (p (ip o T n ) dm — J (p dm J ip dm 

And if ip is constant on local stable leaves then 

cp (ip o T n ) dm — (p dm / ip dm 


<ce n m„ mi , 


Ml. 


( 2 , 2 ) 

(2,3) 


Among C 2 Anosov diffeomorphisms the subset of the ones that admit no invariant measure 
absolutely continuous with respect to m is open and dense. However, it is known that 
any Anosov diffeomorphism of the torus is conjugate to a linear model; see [2S]. The 
conjugacy H (and its inverse) is Holder continuous [21] and carries over the measure m 
to a probability measure // = H*pm) on T 2 which is invariant and mixing by the original 
Anosov diffeomorphism. In particular, the exponential decay of correlations (2.2) is still 
valid for /i, although corresponding to different values of C, 9 and a. We remark that the 
conditions we will analyze in the next section, to ensure the existence of an Extreme Value 
Law, are also invariant under conjugacy. 


The starting point of our analysis is a stationary stochastic process Ao, Ah, Ad,... arising 
from the system described above in the following way. We fix some point ( G T 2 and let 
p : T 2 —* M U {+ 00 } be given by 


p( z ) = - \og(d(z,Q), 

where d denotes the projection of the usual Euclidean metric in 
define X 0 = p and for all i G N 

Xi = tpo T\ 


(2,4) 


to T, as in ( |2.1[ ). We 
(2.5) 


Since m is an invariant probability measure on T 2 , then the stochastic process A"o, Ad, Ad,... 
just defined is stationary. 

Note that {A 0 > u} = B e -n(Q, where B e -u(() denotes the ball centered at ( of radius e~ u , 
in the Euclidean metric. In other words, having an exceedance of an high threshold u, at 
time j, means that, at the jth iterate, the orbit enters the ball around ( of radius e~ u . 

We remark that we could use a function different from — log in the definition of p. In fact, 
any not too irregular function, achieving a global maximum at 0 and invertible in a vicinity 
of 0, would work. See [9j Section 1.1] for conditions on this function. For definiteness, here 
we will use — log which puts us on the Type 1 or Gumbel domain of attraction for the 
maxima. Our goal is to show the existence of an EVL and for that purpose the analysis 
does not change much from one example to another. 
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In order to obtain a non-degenerate limit for the distribution of M n := max{A" 0 ,..., X n _i}, 
one usually uses normalizing sequences (w n ) ne N such that the average number of ex¬ 
ceedances among the first n observations is some constant r > 0, when n is large enough, 
i.e., the sequence ( u n ) n& ^ satisfies: 

nm(X o > u n ) -» r > 0. (2.6) 

n—>■ oo 


In the case here, if we take u n = | log(7rn) — | logr, then conditions (2.6) is satisfied. 

We are now in condition of stating our first result which asserts a dichotomy regarding 
the possible limit distributions for M n depending on the periodicity or aperiodicity of the 
point (. 

Theorem 2.2. Let T : T 1 2 —>■ T 2 be an Anosov linear diffeomorphism. Fix ( e T 2 and 
define the stochastic process X 0 ,Xi,X 2 , ... as in (2.5). Let ( u n ) ne n be a sequence satisfying 
( 2.6[ ) for some r > 0. Then 

(1) if ( is not periodic, we have lim„_> 0 0 m(M n < u n ) = e _r . 

(2) if ( is periodic of prime period q, we have lirn^oo m(M n < u n ) = e~^ T , where 

x) — - ( arcsin , — arcsin , 1 * 

w V Vl A l 29+1 


In contrast to the examples of uniformly expanding systems (with no contracting direc¬ 
tions), the metric used plays an important role in the computation of the EL To illustrate 
this, let e u ,e s denote unit eigenvectors associated to the eigenvalues A and 1/A, respec¬ 
tively; thus {e M ,e s } forms a normal basis of M 2 * . Let 2 = ( x u ,x s ), where ( x u ,x s ) stands for 
the coordinates of z in this basis. Define the metric d* by 

gT(z, 0) = max{|a: u |, |x s |} (2.7) 

and then project it to T 2 * . 


Theorem 2.3. Let T : T 2 —> T 2 be an Anosov linear diffeomorphism. Fix ( G T 2 and 
define the stochastic process Xq,X\,X 2 , ... as in (2.5) for the observable ip as in (2.4) but 
evaluated with the metric d* instead. Let ( u n ) n be a sequence satisfying (2.6) for some 
t > 0. Then 


(1) if ( is not periodic, we have hin^^oc rri(M n < u n ) = e T . 

(2) if £ is periodic of prime period q, we have lim n ^. 00 m(M n < u n ) = e~'° T , where 

d — 1 — | A | — 9 . 


3. The approach and dependence conditions 

In order to prove the existence of EVLs in a dynamical systems context, there are a couple 

of conditions on the dependence structure of the stochastic process that if verified allow 

us to obtain such distributional limits. These conditions are motivated by the conditions 

D{u n ) and D'(u n ) of Leadbetter but were adapted to the dynamical setting and further 
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developed both in the absence of clustering, such as in EQIS1U2], and in the presence of 
clustering in Pi- Very recently, in [12], the authors provided some more general conditions, 
called J\{u n ) and jJ! q (u n ), which subsumed the previous ones and allowed them to address 
both the presence (q > 1) and the absence (q = 0) of clustering. To distinguish these 
conditions the authors used a Cyrillic D to denote them. We recall these conditions here. 


Given a sequence (u n ) nG n of real numbers satisfying (2.6) and q E No, set 

A^"> : = {X 0 > u n ,X 1 <u n ,...,X q < u n }. 

For s, £ G N and an event B, let 

s+l-l 

wsAb)= n t ~\b c )- 


Condition (JX q (u n )). We say that J\{u n ) holds for the sequence X Q ,Xi,. 

£ t Tl G 

where y(g, n, t) is decreasing in t and there exists a sequence (t n )neN such that t n = o(n 
and ny(g, n, t n ) —>• 0 when n -E oo. 


(3.1) 
if, for every 

(3.2) 


For some fixed q G No, consider the sequence (f n )neN given by condition J\ q (u n ) and let 
(k n ) n m be another sequence of integers such that 

k n —y oo and k n t n = o(n). (3.3) 

Condition (JX q ( u n))- We say that JX q {u n ) holds for the sequence X 0 , Xi ,... if there exists 
a sequence (fc n )neN satisfying (3.3) and such that 


|n/fc„J 

lim n P i A n ] n T~ j (A^)) = 0. 

3 = 1 

We note that, when q — 0, condition JJ! q (u n ) corresponds to condition D\u r , 
Now let 

t? = lim d n = lim P ^ 4 wl. 

n->oo n—»oo P(f/ n ) 


(3.4) 


from 


(3.5) 


From m Corollary 2.4], it follows that if the stochastic process X 0 , W,... satisfies both 
conditions and JX q (u n ), the limit in (3.5) exists and 

lim P(M n < u n ) = e^ 0r . 


Hence, the strategy to prove Theorems 2.2 and 2.3 is to show that conditions Ao(wn) and 
JXo(u n ) hold when ( is a non periodic point, and that conditions J\ q {u n ) and JX q (u n ) hold 
when C is a periodic point of prime period q. In fact, we only check these conditions for 
the usual Euclidean metric as in the context of Theorem T2, which is technically harder, 
leaving the necessary adjustments when dealing with the adapted metric d* for the reader. 
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4. Checking condition J\ q {u n ) 


Let flbea set homeomorphic to a ball whose boundary is piecewise smooth, so m(dD) = 0, 
and define 

H k (D) = {x G D : T k (W°(x )) D dD ± 0} . 

Proposition 4.1. There exist constants C > 0 and 0 < 77 < 1 such that, for all k, 

m(H k (D )) < Cr k . (4.1) 

Proof. As a consequence of the uniform contraction within local stable manifolds, there 
exists Ci > 0 such that dist (T n (x),T n (y)) < Ci/|A| n for all y G W*(x). In particular, this 
implies that \T k (Wl(x))\ < Ci/|A| fc . Therefore, for every x G H k (D), the leaf T k (W((x)) 
lies in an annulus of width 2/|A| fc around dD. By Remark 2.1 and the invariance of m, the 
result follows with T\ = 1/|A|. □ 


Lemma 4.2. Suppose <L : M —>■ 


is a Lipschitz map and T is the indicator function 
T := 1 


n,e(A ( n q) )' 


Then for all j > 0 

' <L (T o T- 7 ) dm 


&dm / Tdm 


<C (||<L|| 00 Tp J + ||4> 


I Lip 


Qti/ 2J 


(4.2) 


Proof. We choose for reference the fixed point Q of T and its local unstable manifold 
7 “ := Wf((). By the hyperbolic product structure, each local stable manifold W({x) 
intersects 7 “ in a unique point x. Therefore, for every map T we define the function 
T(a;) := v L(a;). Also, for every *gN, let 47 = T o T l . We note that 4q is constant along 
stable manifolds and that the set of points where v Iq ^ \Eq is, by definition, the set of x for 
which there exist r G No and 27 , £2 on the same local stable manifold as T r (x ) such that 

X! G W i4 (4 9) ) 

but 

^ i (4 9) ) • 

Moreover, this set is contained in 1 H k ( ArP). If we let i > |_j /2J then by Proposi¬ 

tion |4J] we have 

OO 

m({T^T ?: })< Y, rn(H k (A^))<Crl j/2i . 

fc=Li/2j 


By (|2.3|), we also get 

^ <I> ( 4> y 3 y 2 j o T- 7- Lt/ 2 J) dm — J $dm J vj/ ^ ^ 2 j dm 

Using the identity 

j 0 ( , 0 oT)— f <f j if = j cp('ifoT — 'ijjoT)+ j <f{ipoT) 


< C ||4>|| L i P Halloo 0 Vjm - 
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we obtain 

$ ((^ Lj/ 2 j-¥ b/ 2 j)o^-^ 2 J) dm 

dm 


< 


+ C ll$l| Lip ^/ 2J 


$ dm J L// 2J - ^Li/ 2 j) 

< C (2 Halloo m {% /2J ^ Vb u/2J } + ||$|| Lip & W J) 

< c (||$||oo rP J + ||$|| Lip ^ /2J ) . (4.3) 

We complete the proof by observing that J T dm = j 2 j dm, due to the T-invariance 
of m and the fact that v h^y 2 j o T't-U/2j = vjn = VJ/ 0 ju. □ 

To prove condition ,U,(u n ), we will approximate the characteristic function of the set A$ 
by a suitable Lipschitz function. However, this Lipschitz function will decrease sharply to 


zero near the boundary of the set An*. As the estimate in Lemma T2 involves the Lipschitz 
norm, we need to bound its increase as we approach dAn*. 

Let A n = A^ and D n := 16 An' 1 : dist (x, Afy > n~ 2 \ , where denotes the closure 


(4.4) 


of the complement of the set A n . Dehne : X —> M. as 


0 


if x d A r 


* — l dist L'W) If c/1 \ n 

di S t(a;,A®)+dist(x,D n ) 11 x ^-n \ 

1 if x G D n 


Note that <3> n is Lipschitz continuous with Lipschitz constant given by n 2 . Moreover, as 

we have ||$ n — IaJIlHto) A C/n 2 for some constant C. 


observed in Remark 
It follows that 


2.1 


1 A ( q ) (^[j/ 2 j oR dm — m(A$) / Tdm 


< 


1 A (g) - 'b|j/2j dm + C ^||$n||oo f T } j/4J + ||$n||Lip 0 [ 


1 ^( 9 ) - ) dm I T yy 2 j dm 


(4.5) 


and consequently 


|m (Af n ~ MAi?') m (V„ .,(A<?>)) | < 7 (n,j) 

7 (n,j) — C ( J1 ~ 2 + $L^) 


where 
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and 


9\ = max {ti, 9} . 

Thus if, for instance, j — t n — (logn) 5 , then 717 ( 71 , t n ) —> 0 as n —> oo. Note that we have 
considerable freedom of choice of j in order to ensure that the previous limit is zero; taking 
into account the possible applications, we chose t n = (logn) 5 . 


5. Checking condition JXJu n ) 


Before checking the condition JX' q (u n ), we observe that we only need to consider the sum 
in (3.4) up to t n = (log n ) 5 since, by the exponential decay of correlations, the sum of the 


remaining terms goes to 0. More precisely, we use decay of correlations to show that 

|n/fc„J 

lim n Y" m^nr^) = 0. 

n—Kx) * ^ 

j=log 5 n 

Let be a suitable Lipschitz approximation of 1 (,), defined as in 


(5.1) 


4.4 


Then 


l 4 w(l 4 (?) o T 3 )dm 


1 , 0 ) dm 

Sin. 


< 


o T J )dm- i / <h n dm 


dm 


1 , 0 ) dm 


1 (,)(1 ( ,)oT J )clm- / °T J ) dm 


If (logn ) 5 < j < \n/k n J, then, due to the decay of correlations, the first term is upper 
bounded by 

2 


<L n (<L n o T 3 )dm 


4> n dm 


< Cn 2 9 j < Cn 


-2 


if n is sufficiently large. Recalling Remark \2.1\ for the second term we obtain, for n big 
enough, 

2 / « \ 2 


drnj - \^J dm 
Similarly, we estimate the third term as follows 


< Cm(A^\D n ) < Cn~ 2 . 


$ n ( 4 ) n o T J )dm- / 1 (,)(! <,) oR)dm 


< 2 m(A®\D n ) < C n~ 2 . 


Hence equation (5.1) is satisfied. 
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5.1. The periodic case: checking JX q {u n ) for q > 0. Suppose £ is a periodic point 
of minimal period q for the map T. If G = T q , then £ is a fixed point for the Anosov 
diffeomorphism G, which has eigenvalues A anc [ j\9 ; with < 1 < |A |, where j and A 
are the eigenvalues of the original map T. 


la 


Note that An corresponds to the shaded region in Figure 
{X 0 > u n }, while the ellipse is the pre-image T~ q ({X 0 > u n }' 


of An \ we just need to find the distance between the points A and B on the intersection 
of the ball with the ellipse. If we use local coordinates x, y for the unstable and stable 
local manifolds through £, respectively, then the ball’s equation is x 2 + y 2 = s 2 and the 
ellipse’s equation is \\\ 2q x 2 + |A| ~ 2q y 2 = s 2 , where ~ r/(7rn). Then the x coordinate 


the ball represents the set 
To compute the diameter 


of the points A and B is 


yiAI^+i’ 


while the y coordinates are 


|A|<3 


\J\\\ 2q +l 

counterpart, respectively. This means that the diameter of A^ is 


and its symmetric 


y/\X\ 2q +l ' 


Note that, under G~ q , the local stable manifolds in a small ball around £ expand uniformly 
and the local unstable manifolds contract uniformly. By construction of the sets An it 
follows that G~ q (An^nAn^ = 0, while the diameter of G~ q (An^), which is measured along 
the stable direction, is less than 1. Hence, G~i(A$) D A^ = 0 for all j = 1,... ,g(n), 
where 


logn + log(|A| 29 + 1) — 2 log(2 \\\ q (r/Tr) 1 / 2 ) 
2 q log |A| 


Since £ is periodic of period q. for large n we have T q An' 1 D An^ = 0 for j = 1,..., qg(n). 
Thus, 


9 g(n) 

^m(4”nr-J(4’>)) = o, 


(5,2) 
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Let x G ArP Cland 7 x := W^P\An\ For j > qg(n) , the curve T 3 p) x ) may be sufficiently 
long to wrap around the torus. Yet, if this is the case then 

sjtl 

Indeed, the length of APP in the stable direction is bounded above by C\ - 7 =, for some 
constant C\ independent of the local stable manifold. As T restricted to local stable 
manifolds has bounded distortion (in fact its derivative at any point is even constant, 
equal to L ), it follows that \{y : y G 7 ^, T~ 3 {y) G An'*}\ < ^\lx\- Moreover, since the 
2-dimensional Lebesgue measure m is invariant by T, we have 




= m(AMnTi(AM)) 



\ lx nP(A^)\dm x 


c 


c 


\lx\dm x < —;=■ m(A^). 


'Ai q) nw: 


n 


This implies that 

log 5 n 

n n T~ 3 (A^)) = O (n(logn ) 5 n~ 1 ^ 2 m(A^)) = 0 (n _1//2 log 5 n) 

j=qg{n)+l 


which, together with (5.1) and (5.2), yields (3.4). 


-> 0 

n—)-oo 


(5.3) 


5.2. The non-periodic case: checking JX q {u n ) for q = 0. Let us consider now a non¬ 
periodic point (. Note that, in this case, An'* = U n = {A " 0 > u n } is a ball centered at ( 
such that m(U n ) ~ r/n. We split the sum in (3.4) into three parts: 


\n/k n \ g(n ) log 5 n 

n m{U n C\T-i(U n )) = n^m{U n nT-i(U n )) + n m(U n n T~ 3 (U n )) 


3 =0 


3 =0 


j=g(n )+1 


n/kn 


+ n Y, m(U n nT- j (U n )), 


(5.4) 


j=log 5 n +1 


where g(n) = 1S such that, for j > g{n ), the set T 3 (U n ) is sufficiently stretched 

along the stable direction to start wrapping around the torus. The third summand on the 
right of (5.4) goes to 0 on account of ( |5.1 ). We proceed considering the other summands, 
starting by the first. Contrary to the periodic setting, we now have to deal with the 
possibility of U n D T~ 3 (U n ) 7 ^ 0 for j < g(n). Let R n = min{j G N : U n D T~ 3 {U n ) 7 ^ 0}. 
By continuity of T and the fact that A n U n = {£}, we have that R n -A 00 as n —)■ 00 . 


For j < g(n), observe that T 3 (U n ) is an ellipse centered at T 3 (C), stretched along the 
stable direction by a factor Ap and contracted on the unstable direction by l/|Ap. Also, 
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notice that, as we are dealing with a projection on the torus of a linear map in M 2 , the 
stable and unstable directions are always lined up regardless of the position of 
Now, since the set T~ 3 (U n ) is not long enough (on the stable direction) to wrap around the 
torus, the diameter of U n is 2 s n , where s n ~ y/r / y/wn. The length of T~ 3 {U n ) along the 
stable direction (which coincides with the diameter of T~ 3 {U n )) is |Ap's n ; and the length 
of T~ 3 {U n ) on the unstable direction is | A | s n . So, by the geometry of the sets, the area 
of the intersection U n D T~ 3 (U n ) occupies a portion of T~ 3 {U n ) corresponding at most to 
a constant times |A| _J of its area. That is, there exists some C > 0 independent of j such 
that 

m(U n n T-*(U n )) < C |A|-J m(T-i(U n )) = C |A|^‘ m(U n ). 

In the case of the Euclidean metric we can take C ~ 2/ 7 T. Then it follows that 

g(n) g(n) 

nV m(U n DT- j (U n )) =n V \X\- j r/n < Ct |A|- r " -> 0. 

‘ ‘ n—toa 

j =0 j=Rn 


We are left to estimate the second summand of (5.4). Let x G U n nW^ and 7 ^ := W* fl U n . 
For j > g(n), the curve T~ 3 (j x ) may be sufficiently long to wrap around the torus. Set 


Ui ■■= U 7x- 

{x e u n nw™ : |t-j(7x)I>i/2} 

Note that for the usual Euclidean metric we have 1E“ fl U n \ U 3 n \ < 2a/s 2 — |A|~ 2 t/16. 
Since for x G U 3 V we have \T~ 3 ('y x )\ >1/2 and the diameter of U n is less than 2 s n , then 
there exists C > 0 such that, for all j > g(n ) and all 1 G f/j, we have 

\T~ j ^ x )DU n \ <Cs n \T~ 3 { lx )\. 


As T restricted to local stable manifolds has bounded distortion (since the derivative is 
constant), it follows that \{y : y G 7 ^, T~ 3 (y) G U n }\ < Cs n | 7 X |. Moreover, by the 
T-invariance of m we have 

m{U n n T-*(U£)) = m(Ui n T 3 (U n )) = ! 1 7 a: n T 3 (A^)\dm x 

Juinw“ 

< C s n f 1 7a: I dm x <C s n m{U 3 n ). (5.5) 

■Julnw™ 

Now, we focus on T~ 3 (U n \U^) which is a made of two connected components with diameter 
1/2 along the stable direction and with length |A| -3 s n ^l — \Jl — |A| _2 l/(16s ri 2 )j. As 

before, the geometry dictates that the area of the intersection U n fl T~ 3 (U n \ 17/) occupies 
a portion of T~ 3 {U n \ 17/) that corresponds at most to a constant times s n of its area. 
Indeed, there exists C > 0 such that, 

m(U n n T-l(U n \ Ui)) < Cs n m(T~ 3 (U n \ C7/)) = C s n m(U n \ U 3 n ). (5.6) 
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Joining (5.5), (5.6) and recalling that s n ~ y/r/y/im, we deduce that 

log^n 

n y m{U n fl T~ j (U n )) = 0{n{\ogn) b n~ l ^ 2 m{U n )) = 0(n -1 / 2 log 5 n) -* 0. (5.7) 

j=g(n)+l 


6. Computing the Extremal Index 


We recall that the Extremal Index $ is given by the limit in (3.5). 


In the case of a non-periodic point £, we have q = 0 and the computation is trivial because 
An' 1 = U n , hence d — l. 


If ( is a periodic point of prime period q, we need to compute the area of the shaded 
regions in Figure la Recalling that the second coordinate of the point A of the picture is 
y A = -^== i and using the symmetry of the region, we have 


fVA 


m Oi !, ) = 4 / V s S-r ! Ph 2 ' , % 


( 6 . 1 ) 


(. 


^ n arcsin 




Thus, in this case, 


= — arcsin 


f. 


\/\\\ 2q + 1 
|A|« 


— arcsm 


\/\X\ 2q + 1 


7r 




\/\X\ 2q + 1 


arcsm 


viwir 


( 6 . 2 ) 


Remark 6.1. It is interesting to observe that the metric used affects the value of the 
Extremal Index. In fact, if we had used the adapted metric, so that its balls would 
correspond to squares with sides lined up with the stable and unstable directions, then 
the El would be equal to 1 — |A | ~ <? , which is unsurprisingly more in tune with the one- 
dimensional setting. See Figure [lb] Notice also that, in both metrics, the extremal index 
approaches the non-periodic value as the period goes to oo, that is, lim^+oo = 1. 


7. Conclusion 

For C 2 expanding maps of the interval a dichotomy exists between the extremal behavior of 
periodic points and non-periodic points 0 mm- In this paper we demonstrated that the 
same dichotomy holds for hyperbolic total automorphisms. The extreme value statistics of 
observations maximized at generic points in a variety of non-uniformly hyperbolic systems 
is known to be the same as that of a sequence of iicl random variables with the same 
distribution function m flOl I9l H6j but the statistics of functions maximized at periodic 
points has not yet been established. We believe that it will be the same as in the iid case 
for smooth non-uniformly hyperbolic systems but expect mixed distributions (as in HI) for 
non-uniformly systems with discontinuities. 
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8. Rare Events Point Processes 


In order to make the presentation more tractable we started by stating the dichotomy in 
terms of the possible limit distributions for the maxima when the centers are chosen to be 
either non-periodic or periodic points. However, we can enrich the results without much 
extra work by considering the convergence of point processes of exceedances or rare events. 

If we consider multiple exceedances we are led to point processes of rare events counting 
the number of exceedances in a certain time frame. For every icKwe define 

A^(A) := ^2 lxi>u- 

ieAnNo 

In the particular case where A = [a, b) we simply write A^ a := A<u([a,b)). Observe that 
counts the number of exceedances amongst the first n observations of the process 
X 0 ,X\,... ,X n or, in other words, the number of entrances in U{u ) up to time n. Also, 
note that 

{•'Co = 0} = {M„ < u} (8.1) 


In order to define a point process that captures the essence of an EVL and HTS through 
(8.1), we need to re-scale time using the factor v : = 1 /P (X > u ) given by Kac’s Theorem. 
Let S denote the semi-ring of subsets of R^ whose elements are intervals of the type [a, b), 
for a, b G R(j~. Denote by IZ the ring generated by S ; recall that for every J G IZ there 
are k G N and k intervals I{, . . . , Ik £ <S , say Ij = [a 3 , b 3 ) with aj,bj e Rq, such that 
J = UFor I = [a, b) € S and a G R, we set al := [aa,ab) and I + a := [a + a,b + a ). 
Similarly, for J G IZ, we define aJ := al\ U • • • U ah and J + a := (h + a) U ■ ■ ■ U (h + <r). 


Definition 1. Given J ElZ and sequences (u n )n.eN and (u n )neN, the rare event point process 
(REPP) is defined by counting the number of exceedances (or hits to U{u n )) during the 
(re-scaled) time period v n J G IZ. More precisely, for every J ElZ, set 


N n (J) A<a n (v n J) — ^2 lXj>u n - 

j€v n JnNo 


( 8 . 2 ) 


Under dependence conditions similar to the ones previously analyzed, the REPP just de¬ 
fined converges in distribution to a standard Poisson process when no clustering is involved, 
and converges in distribution to a compound Poisson process with intensity d and a geo¬ 
metric multiplicity distribution function otherwise. For the sake of completeness, we now 
define what we mean by a Poisson and a compound Poisson process. (See |20j for more 
details.) 

Definition 2. Let Yf, I 2 , ■ • • be an iid sequence of random variables with common exponen¬ 
tial distribution of mean l/th Let Z 1 , Z 2 ,... be another iid sequence of random variables, 
independent of the previous one, and with distribution function n. Given these sequences, 
for J <ElZ, set 

1 j d 
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where 5 t denotes the Dirac measure at t > 0. Whenever we are in this setting, we say that 
N is a compound Poisson process of intensity d and multiplicity d.f. n. 

Remark 8.1. In this paper, the multiplicity will always be integer valued. This means that 
7r is completely defined by the values 7r fc = P {Z\ = k ), for every k G N 0 . Note that, if 
7Ti = 1 and -d = 1, then N is the standard Poisson process and, for every t > 0, the random 
variable N([0,t)) has a Poisson distribution of mean t. 

Remark 8.2. When clustering is involved, typically 7r is a geometric distribution of param¬ 
eter •& G (0,1], i.e., TTk = $(1 — d)^ 1 , for every k G No- This means that, as in [IT], the 
random variable iV([0,t)) follows a Polya-Aeppli distribution 

P(JV([0, *))) = * = e-* ^(1 - ( k I 1) ' 

j =1 J' \J J 

for all fcGN and P(JV([0,t))) = 0 = e - **. 


We are now in condition to state the dichotomy in terms of the convergence of REPP. 

Theorem 8.3. Let T : T 2 — y T 2 be an Anosov linear diffeomorphism. Fix ( G T 2 and 
define the stochastic process A 0 , W, X 2 ,... as in ( 2.5| ). Let (w n )neN be a sequence satisfying 
( 2.6[ ) for some r > 0 and (u n )neN be given by v n = l/P(Ao > u n ). Consider the REPP N n 
as in Definition [2j Then 

(1) if f is not periodic, the REPP N n converges in distribution to the standard Poisson 
process. 

(2) if f is periodic of prime period q, the REPP N n converges in distribution to a 
compound Poisson process with intensity d given by (6.2) and multiplicity d.f. n 
given by (8.6). 


As in the computation of the El, the metric used has an important impact on the multiplic¬ 
ity distribution obtained when we have convergence of the REPP to a compound Poisson 


process. Indeed, a similar result for the adapted metric d* given in (2.7) is as follows 


Theorem 8.4. Let T : T 2 —>■ T 2 be an Anosov linear diffeomorphism. Fix ( G T 2 and 
define the stochastic process X 0 , X 1 , X 2 ,... as in (2.5) for the observable in ( |2.4[ ) but 
evaluated with the metric d* instead. Let (u n ) n ^ be a sequence satisfying ( |2.6[ ) for some r > 
0 and (ivJngN be given by v n = 1/P(A 7 " 0 > u n ). Consider the REPP N n as in Definition ^ IJ 
Then 


(1) if f is not periodic, the REPP N n converges in distribution to the standard Poisson 
process. 

(2) if f is periodic of prime period q, the REPP N n converges in distribution to a 

compound Poisson process with intensity d = 1 — |A| _<? and geometric multiplicity 

d.f. tv* given by i t*(k) = 19(1 — i?)” -1 , for all kGN. 
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8 . 1 . Absence of clustering. When condition JX Q {u n ) holds, there is no clustering and so 
we may benefit from a criterion, proposed by Kallenberg [ 2 D, Theorem 4 . 7 ], which applies 
only to simple point processes without multiple events. Accordingly, we can merely adjust 
condition ,H,o(w n ) to this scenario of multiple exceedances in order to prove that the REPP 
converges in distribution to a standard Poisson process. We denote this adapted condition 
by: 

Condition (. D 3 [u n )). Let A E 1 Z and t E N. We say that D 3 (u n ) holds for the sequence 
Xo, X ± , ... if 

|P ({Xo > u n } n {jV Un (A + t) = 0}) - P({X 0 > u n })F(^ Un (A) = 0)| < 7 (n, t), 

where 7(71, t) is nonincreasing in t for each n and 717(71, t n ) -)■ 0 as n -)• 00 for some 
sequence t n = o(n). (The last equality means that t n /n —* 0 as 71 —* 00). 


In P Theorem 5] it was proved a strengthening of Theorem 1] which essentially says 
that, under conditions D 3 (u n ) and JX 0 (u n ), the REPP N n defined in (8.2) converges in 
distribution to a standard Poisson process. 


The proof of condition D 3 (u n ) follows after minor adjustments to the proof of JX q (u n ) in 
Section]!] Since condition jX a {u n ) holds at every non-periodic point ( (see Section [5(2]), then 
for all such points £ the corresponding REPP N n converges in distribution to a standard 
Poisson process. 


8.2. Presence of clustering. Condition jX Q {u n ) prevents the existence of clusters of 
exceedances, which implies that the EVL is standard exponential H(t) = e~ T . When 
JX q {u v ) fails, the clustering of exceedances is responsible for the appearance of a parameter 
0 < 7? < 1 in the EVL, called El, and implies that, in this case, H(t) = e -l9r . In [TO], the 
authors established a connection between the existence of an El less than 1 and a periodic 
behavior. This was later generalized for REPP in [ITT) . 

For the convergence of the REPP when there is clustering, one cannot use the aforemen¬ 
tioned criterion of Kallenberg because the point processes are not simple anymore and 
possess multiple events. This means that a much deeper analysis must be done in order to 
obtain convergence of the REPP. This was carried out in m and we describe below the 
main results and conditions needed. 

Let C be a periodic point of prime period q. Firstly, we consider the sequence (£/( K - ) ('U ri )) K>0 
of nested balls centred at ( given by 

l/ w) (ti„) = U(Un) = U„ and [/<“>(«„) = n U(u n ), for all k £ N. 

(8.3) 

Then, for i, k, £, s E N U {0}, we define the following events: 

:= T-‘ ([/<">(«„) - t/<“ +1) (ti„)) . (8.4) 

Note that Q° 0 (u n ) — AnK Besides, U n — IJ)1 0 Q * 0 (u n ), which means that the ball centred 
at ( which corresponds to U n can be decomposed into a sequence of disjoint strips where 
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FIGURE 2. The sets Q q0 {u n ) for the Euclidean metric. 


Qq,o( u n ) are the most outward strips and the inner strips Q q ^ l {u n ) are sent outward by T p 
onto the strips Q qfi (u n ), i.e., T q (Q q ^ 1 (u n )) = Q qfl (u n ). 

We are now ready to state the adapted condition: 

Condition (. D q {u n )*). We say that D q (u n )* holds for the sequence X 0 , Xi, X 2 ,.. ■ if for 
any integers t, fiq,..., k^, n and any J = U j =2 Ij £ H with inf{x : x G J} > t, 

|Pn (n= ^)) - p (<?"„(«„)) p = Kj )\ < 7 («, J). 

where for each n we have that 7(n, f) is nonincreasing in t, and 717(71, —> 0 as n —y 00, 

for some sequence t n = o(n). 


In m, for technical reasons only, the authors introduced a slight modification to JX' q (u n ). 
This condition was denoted by D' q {u n )* and it requires that 

[n/k n ] 

lim n V P(Q, 0 (wn) H {X, > u n }) = 0, 

n—> 00 z 

i=1 


which holds whenever condition does. 

From the study developed in it follows that if X 0 , X \,... satisfies conditions D q (u n )*, 
D' q {u n )* and lim^oo X) K >i P(C (k) (u n )) = 0, where ( u n ) nG ^ is such that ( |2.6[ ) holds, then 
the REPP N n converges in distribution to a compound Poisson process with intensity d 
given by (3.5) and multiplicity d.f. 7 r given by: 

(p«;.o 1 k))-pw;,o(«o)) 


7 t(k) = lim 


P(Q%M) 


(8.5) 


8.3. Computing the multiplicity distribution. In order to compute the multiplicity 

distribution, we need to calculate m(Q q0 ). To do so, let E % denote one of the four inter¬ 
section points of dT~ iq (U n ) with dU n . Let t/eJ denote the coordinates of E t on the 
unstable and stable directions, which correspond to the vertical and horizontal directions, 
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w“ 


w s 


Figure 3. The sets Qg 0 (u n ) in the adapted metric. 


respectively, in Figure [2} We assume that the choice of the Ei s among the possible four 
points is made so that xE t , yE. t > 0. Hence, 

Sn|A| i<? 


Ve z - 


\/T + X 2iq ’ 

because dT~ iq (U n ) is described by the equation X 2iq x 2 + X~ 2iq y 2 = s\. Therefore, for n > 1, 
we may write 


"VEi 


m 


(Qq,o( u n)) = 4 / A - X~ 2iq y 2 dy + 4 / yjsi - y 2 dy 


yE i+ 1 

VE f 


niE. 


i +1 


-4 / S 2 n - X-W+I) q y 2 dy 

|A|h+ 1 )9 . 1 


= 2 si ( arcsin ._ 

V y/l + X 2( x +V > q 


arcsm 


+ 2 si ( arcsin . = 

n V y/l + A 2 *9 


arcsm 


y/l + A 2 (* +1 )9 

|A|* 9 


y/l + X 2iq , 

Recall that Qqo(u n ) = An^ and its measure has been computed in (6.1). Hence, using 
formula (8.5), we obtain 


7TIK = 


( ai vom / - 

a/ i+a 2 ^- 1 ^ 


arcsm 


— arcsm 


lAlU-m \ , ( ■ 1 ■ | AI C«+i)« 

—- 1 + arcsm-° ,T1 ——- 


\J l+A 2 ('=- 1 )9 J Y a/ 1 +A 2 ( k +!) 


— arcsm 


A/ 1 +A 2 ( k + 1)9 


arcsin tot - arcsin TOT 


| 2 (arcsin arcsin 

arcsin TOT “ arcsin TOT 


( 8 . 6 ) 


Remark 8.5. As for the EI (see Remark 6.1), the multiplicity distribution depends on the 
metric used. If we were to use the adapted metric d*, in which balls correspond to squares 
whose sides are lined up with the stable and unstable directions, then we would obtain a 
geometric distribution for the multiplicity, i.e., 7t*(k) = T*(l — T*) K_1 , where d* = 1 — |A|“ <3r . 
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Observe also that in the Euclidean metric case we do not obtain a geometric distribution, 
as one did for repelling periodic points (see m or in the case of the adapted metric 

However, taking into account that 


mentioned in Remark 8.5 


arcsm 


x 


vTT 


x* 


7T 1 

—- Vo{l/x) and arcsin 


x 


vTT 


x i 


= - + o( 1/x), 

X 


one deduces that 


ii m hYth = iA|-’ = ,r,(K + 1 ) 


K^-oo 7 t(a) 


7 T*(k) 


where is as given in Theorem 8.4 Hence, we obtained a multiplicity distribution 

different from the geometric distribution but which is not that far from it. 

The proof of condition D q {u n )* follows with some adjustments from the proof of J\ q (u n ) 
in Section 4| Since condition JX q {u n ) holds at every periodic point £, as was shown in 
Section|5.2 then condition D' q (U n )* from |Tl] also holds for all such points £. Let C% denote 
a rectangle, centered at £, with its sides lined up with the stable and unstable direction, 
of length s n on the stable direction and |A|“ KIJ s n on the unstable direction. Observe that 
UM (u n ) C C£, which implies that 


m(U^\u n )) < |Ar K< 0 


and consequently 


lim Vm(£/ ( 

n .—yon * ^ 


k) 


Ur 


K> 1 


= lim 

n—¥ oo 


:st = 0 . 


1 “ |A|-« 

So, at every periodic point ( of prime period q, the respective REPP N n converges in 


distribution to a compound Poisson process with intensity t? given by (6.2) and multiplicity 


d.f. 7T given by (8.6). 
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